Scalar spin chirality and quantum Hall effect on triangular lattice 
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We study the Kondo Lattice and Hubbard models on a triangular lattice for band filling factor 3/4. 
We show that a simple non-coplanar chiral spin ordering (scalar spin chirality) is naturally realized 
in both models due to perfect nesting of the Fermi surface. The resulting triple-Q magnetic ordering 
is a natural counterpart of the collinear Neel ordering of the half-filled square lattice Hubbard model. 
We show that the obtained chiral phase exhibits a spontaneous quantum Hall-effect with a^y ~ e^//i. 

PACS numbers: 71.10.Fd, 71.27.-f-a, 73.43.-f 
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Interactions commonly lead to magnetism. For in- 
stance, the ground state of the half-filled Hubbard model 
becomes antiferromagnetic (AFM) for a finite onsite 
Coulomb interaction U. The AFM phase is collinear on 
bipartite lattices and it exists for any finite U/t {t be- 
ing the nearest-neighbor hopping) in the cases of square 
and cubic lattices due to nesting properties of the Fermi 
surface (FS). The situation is less clear when compet- 
ing interactions or lattice geometry lead to frustration 
and more complex magnetic orderings. For instance, a 
120° three-sublattice (non-coUinear) structure is stabi- 
lized on the half-filled triangular lattice for large enough 
U/t and a spin liquid ground state has been proposed for 
the kagome lattice. The non-coplanar spin orderings have 
proven to be particularly elusive. Besides being the next 
level in magnetic ordering complexity, the non-coplanar 
states have unusual transport [l[ and magneto-electric 
properties resulting from their chiral structure. Conse- 
quently, it is important to find simple and realistic models 
in which such states could emerge and remain stable over 
a range of Hamiltonian parameters. 

A triangular plaquctte is a basic building block for lat- 
tices with built-in geometrical frustration. Assume that 
on every site i there is an ordered magnetic moment Si. 
If the moments on the plaquette are non-coplanar, the 
resulting scalar spin chirality, (Si • Sj x S^) 7^ 0, breaks 
both the time- reversal and the parity symmetries. No- 
tably, the scalar spin chirality can exist even without 
having statically-ordered magnetic moments ((Si) = 0) 
[1]. When conduction electrons propagate through such 
a spin texture, they may exhibit spontaneous Hall effect 
in the absence of any externally applied magnetic field 
Moreover, if the spin ordering opens a full gap in 
the charge excitation spectrum, a spontaneous quantum 
Hall insulator can form [l| . This exotic phenomenon has 
the following origin: Due to the exchange interaction, 
the underlying local moment texture aligns the spins of 
the conduction electrons inducing a Berry phase in the 
electron wavefunction, which is indistinguishable from a 
"real" magnetic fiux. Ohgushi et al. [l|], considered a 
model of 2D kagome ferromagnet where spin orbit inter- 
action tilts the local moments so as to produce non-zero 
chirality. The spin configuration corresponded to total 



zero fictitious flux through the unit cell, while still cre- 
ating non-zero quantized Hall conductivity at filling 1/3. 
The authors were primarily concerned with the double 
exchange model, which assumes two species of electrons 
~ the localized ones and the itinerant ones, but also sug- 
gested that similar physics may occur in a single band 
Hubbard-like model in the presence of spin-obit interac- 
tion. Appearance of quantum Hall effects in systems with 
spontaneous bond ordering has recently been discussed in 
Raghu et al. 0]. 

In the present work wc consider a model of itinerant 
electrons interacting with a chiral magnetic state on a 
simple triangular lattice. The magnetic ordering has the 
four-site unit cell, with the local spin orientations corre- 
sponding to the normals to the faces of a regular tetrahe- 
dron, see Fig. [1] This magnetic structure has been ear- 
lier proposed as a likely candidate for a low temperature 
phase of Mn monolayer on Cu(lll) surfaces Q and a pos- 
sible nuclear spin ground state of a two-dimensional '^He 
[3 . We demonstrate that in the presence of itinerant elec- 
trons, this structure naturally appears as a weak-coupling 
instability at 3/4 filling due to the nesting properties of 
the electronic Fermi surface, in direct analogy with the 
magnetic instability in the square-lattice Hubbard model 
at half filling. In the insulating state at 3/4 filling the 
model exhibits spontaneous quantum Hall effect. 

Let us first consider the Kondo-Lattice Hamiltonian, 



H 



J 



(1) 



which describes electrons hopping between the nearest- 
neighbor sites of triangular lattice interacting via on-site 
exchange coupling J with the localized moments Si . The 
operator c^^ (c^^,) creates (destroys) an electron with spin 
a on site i and (Tq,^ = (c^^, cr^^, cr^^) is a vector of the 
Pauli matrices. This model may describe various d and 
/-electron compounds, with J being the ferromagnetic 
Hund's coupling in the first case and the antiferromag- 
netic Kondo exchange in the latter. We consider the 
classical limit S ^ 1 for the localized spins. Then, the 
magnitude of |Si| = S* is fixed and the only remaining 
degree of freedom is the spin orientation. The sign of 
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J is irrelevant (the spectrum does not depend on this 
sign and the eigenstates are connected by a global time 
reversal transformation of the localized spins). In the 
following, we take J > 0. 

The simplest limit of H is the strong coupling, JS/t ^ 
1. Two disconnected sectors appear in this regime: the 
low (high) energy subspace of states in which conduction 
electrons hop between neighboring sites remaining paral- 
lel (antiparallel) to the local moment S^. Within each sec- 
tor, the conduction electron spin orientation is enslaved 
to its position, and the effective hopping amplitude ac- 
quires a quantum mechanical phase [8| . This Berry phase 
is indistinguishable from the Aharonov-Bohm phase in 
an externally applied magnetic field 0. The effective 
Berry flux through a closed path is half of the solid angle 
subtended by the electron spin while moving along the 
path ■ The fluxes experienced by electrons with spins 
aligned and anti-aligned with the local moments are equal 
in magnitude but opposite in sign. For the specific case of 
"tetrahedral" ordering. Fig. [H the sohd angle subtended 
by local moments in corners of each elementary triangu- 
lar plaquctte is clearly Q, = 47r/4. Thus, in the limit of 
strong coupling, the problem is equivalent to one of elec- 
trons on a triangular lattice in a uniform spin-dependent 
magnetic field corresponding to flux (p^ = aix jl per tri- 
angular plaquette. The overall energy splitting between 
the spin up and spin down sectors is 2JS. The energy 
spectrum can be easily found. 



ek = ±JS±. (cos^ kai 
v3 



cos^ ka2 + cos^ kaa 



>l/2 



(2) 



where ai and a2 real space lattice vectors and aa = 
ai — Sl2- Since the real space (magnetic) unit cell has 
doubled, this dispersion is defined on halj of the origi- 
nal reciprocal unit cell. The amplitude of the hopping 
is reduced by a factor \fZ due to the angle Q between 
the neighboring spins, t tcos{0/2) = t/^/i. There 
arc total four energy-split bands - (2 spins) x (2 sites 
in magnetic unit cell). Thus for fillings that are integer 
multiples of 1/4, the system is an insulator. Filling 1/2 
corresponds to full filling of the spin-up sector, and there- 
fore it is equivalent to a band insulator, with u^y = 0. 
The "interesting" fillings are 1/4 and 3/4. From the ex- 
plicit treatment of spinless electrons on triangular lattice 
[llj . one finds that \(J.j.y \ = /h in both cases, but with 
opposite signs, since spin up and spin down electrons ex- 
perience opposite effective magnetic fluxes. 

Now we turn to the case of an arbitrary JS/t. From 
Fig. [U the unit cell in the ordered state is quadrupled - 
it is now spanned by vectors 2ai and 2a2. It is conve- 
nient to treat to problem in the momentum space. For 
that we notice that the "tetrahedral" magnetic ordering 
is equivalent to the so-called 3g non-coplanar state, with 
the order parameter 



M =(cos Q|r, cos Q^r, cos Q3r) 
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FIG. 1: (Color online) Chiral spin ordering on a triangular 
lattice, (a) Four possible orientations of the local magnetic 
moments. They correspond to the normals to the sides of a 
regular tetrahedron, (b) Arrangement of the local moments 
on the real space lattice, ai and a2 are the lattice vectors, (c) 
Momentum space. Unshaded hexagons define the Brillouin 
zones. G is one of the reciprocal lattice vectors. Vectors 
Qi are the ordering vectors for the "tetrahedral" chiral order 
{Qi = G/2). Note that Qi perfectly nest the 3/4-fiIled FS 
(shaded hexagons). 



The amplitudes Sa are equal in magnitude but can have 
different signs (arbitrary global rotation of the triad in 
the spin space yields an equivalent state). The ordering 
vectors Qq are equal to half of the reciprocal lattice vec- 
tors, see Fig. [T];. Therefore, momenta -l-Q,, and — Q,j 
are equivalent (Q,, ^ — Qr;), and ±Qa ± Qh ± Qc = 0. 
The scalar chirality, defined as a counter-clock- wise triple 
product Si • (S2 X S3) (X SaSbSc around a plaquette (e.g. 
123 in Fig. [T})) can be either positive or negative, de- 
pending on the choice of Sa- To be specific, let's chose 
all Sa > 0, Fig.[T|i. In momentum space the Hamiltonian 
Eq. ^ becomes 



lea 



-JY. ^vcL'^lpCi.+Q^f^, (4) 



kcK/377 



where r] = a,b, c. The momentum summation is over the 
quarter of the original reciprocal unit cell. The Hamilto- 
nian splits into two disjoint parts 



H = 



H 



c'jHcj + c]jHcii 

f Ck JSa iJSb JSc \ 

Sa £k+Qi --JSc —iJSb 

— iJSb —JSc Ek+Qj JSa 

JSc iJSb JSa Ek+Qa / 



("^kT' "^k+Qii' '^k+Qai' "^k+Qat) 



(5) 
(6) 



). Therefore, the spec- 



S(r) — {Sa cos QaT, Sb cos Qfcr, Sc cos Qcr). (3) trum is (at least) doubly degenerate. The physical reason 
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for the degeneracy is the same as for the degeneracy be- 
tween the spin-up and spin-down band in a mean field 
treatment of an antiferromagnet on a square lattice ~ the 
spin-up and spin-down wave functions are related by a 
lattice translation combined with a spin rotation. Indeed, 
from Fig. [TJ), the local moment arrangement changes as 
(1 <-> 2, 3 <-> 4) under translation by ai. From Fig. [Iji, 
a spin rotation by angle tt around x axis restores the 
initial state. Therefore, the composition of both trans- 
formations leaves H unchanged, while c// ^ c/. 




FIG. 2: (Color) Band structure corresponding to Eq. ([5}. The 
coupling strength is JS/t = 1. The axes are chosen along Qa 
and Qi, and point (0, 0) corresponds to the center of the first 
Brillouin zone. 

Fig. [2] shows the band structure corresponding to H 
in the reduced (quarter) reciprocal unit cell with the 
axes chosen along and Q^. For an arbitrary value 
of JS/t, there are four pairs of degenerate bands. For 
large enough JS/t, there is a gap between all pairs, and 
the dispersion is well described by Eq. The Hall 

conductivity can be calculated directly using the Kubo 
formula. 
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["F(erik) - n-F(ernk)], 



(7) 

with summation performed over all 8 bands n, m and 
ijx,jy) being the x and y components of the current op- 
erator. At zero temperature, when the Fermi distribution 
functions np are 1 for occupied bands and for the empty 
ones, we recover the results for axy obtained earlier for 
JS/t ^ 1. Each pair of degenerate bands contributes 
Chern number ±1, instead of or ±2, due to the addi- 
tional local (in k-space) SU{2) symmetry implied by the 
degeneracy [12 1. 

For small JS/t, the lower three pairs of bands overlap. 
However, the gap between the two upper bands persists 
for arbitrarily small JS/t. This signifies a weak cou- 
pling instability towards 3q ordering at 3/4 filling (1.5 



electrons per site). Indeed, at this filling the FS of non- 
interacting electrons is a regular hexagon inscribed into 
the first Brillouin zone (see Fig. [TJ;). The vectors 
perfectly nest all sides of the FS, and hence a local inter- 
action will cause the 5q instability. In general, an order- 
ing Sa cos(QQ,r) with arbitrarily directed Sa is expected. 
However, the constraint of same magnetization magni- 
tude on every lattice site |Si| = 5* requires that Sa be 
mutually orthogonal, leading to Eq. ([3|). The quantum 
Hall conductivity at 3/4 filling for arbitrarily small JS/t 
is axy = /h at zero temperature. 

We now consider the possibility that the local moments 
do not originate from localized electrons on different or- 
bitals, but are generated sclf-consistently in a single-band 
Hubbard model -ffnubb on a triangular lattice. The con- 
nection between this case and the one studied above is 
evident from the mean-field expression of -ffnubb '■ 



Hubb 



(8) 



Now all the spin degrees of freedom are provided by 
the same (itinerant) electrons and (Si) = (c-^<Tc^Ci^). 
The result of the self-consistent calculation for the non- 
coplanar order parameter as a function of temperature 
is presented in Figure [H The finite-temperature calcu- 
lation of the transverse conductivity based on Eq. ([7|) 
reveals that axy is a good proxy of the order parameter. 
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FIG. 3: Hubbard model, U = 4t. The magnitude of the order 
parameter for non-coplanar 3q structure and Hall conductiv- 
ity as a function of temperature. 

^^Hubb has the same form as Eq. ([T]) but (S,)'^ is no 
longer fixed due to quantum fluctuations (Si is now in 
the S = 1/2 quantum limit). The relaxation of this 
constraint allows for more general 3q structures such 
as coUinear or coplanar. The energy difference between 
the 3q collinear, coplanar, and non-coplanar structures is 
very small, since it is caused by the interaction between 
the order parameters Sa at the corners of the FS. Numer- 
ically we find that the coplanar (but non-collincar) phase 
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]J. The 



is not favored, while the non-coplanar and colKnear [l5| 
3q structures are are very close in energy at the mean field 
level. An analysis beyond mean-field would be required 
to determine the relative stability of the two phases. 

A spontaneous breaking of a continuous symmetry is 
impossible at finite temperature in 2D for short-range in- 
teractions Does this imply that finite-temperature 
spontaneous quantum Hall effect is impossible? To an- 
swer this question we need to analyze the symmetry 
of the problem in more detail. To that end we con- 
sider the Ginzburg-Landau theory for the 3q magnetic 
structure. It can be constructed in close analogy to 
the theory for the Hciscnbcrg antifcrromagnet (AF) on 
a square lattice. However, instead of a single mean- 
field order parameter at momentum (tTjTt), now there 
are three order parameters, with the ordering wavevec- 
tors Qq. To the lowest order, they do not interact 
and each is described by the standard AF Lagrangian, 
Cc. ^ [r + LU^~ (q- Q„)2]|S^ 
propriately scaled frequencies and momenta 
transition is controlled by the parameter r which can be 
tuned by, e.g. temperature. Given the momentum space 
structure of the order parameters, the first symmetry- 
allowed term that selects the mutual orientation of the 
order parameters is Cmt = w[Si • (S2 x 83)]^. Negative 
value of V favors non-coplanar configuration with Sq. or- 
thogonal to each other, while w > favors configuration 
with Si collinear [isj . This term reduces the symmetry 
of the Lagrangian down to 50(3) x Z2; however, being 
higher order in 5, it docs not modify the mean-field tran- 
sition temperature. Non-coplanar ordering corresponds 
to complete breaking of this symmetry. Thus, by the 
number of SO{3) generators, there are three gapless (lin- 
early dispersing) Goldstonc modes. In addition, there are 
three modes that have a small gap ^ —vr^ ju^ induced by 
the orientation-locking term. 

The gapless modes are responsible for the destruction 
of the long-range magnetic ordering at finite tempera- 
ture. However, smooth distortion of the magnetic tex- 
ture does not necessarily destroy the chirality The 
chirality is a discrete (Ising) order parameter, which can 
be destroyed in two ways: Through point defects - vor- 
tices, and line defects - domain walls. Since the order 
parameter space (the "vacuum manifold" ) of the model 
is given by 5*0(3), there are only two topologically dis- 
tinct kinds of vortices - the trivial and the non-trivial 
one 17|, in contrast to the more familiar C/(l) vortices, 
which can have an arbitrary integer winding number. At 
the vortex core, one or more of the order parameters Sq, 
is completely suppressed. The vortices effectively reduce 
the magnitude of average chirality, without changing its 
sign. On the other hand, a domain wall corresponds to 
a line where the sign of chirality changes via mutual re- 
orientation of the order parameters, without necessarily 
changing their magnitude (at the domain wall the state is 



coplanar). The finite temperature transition can be gov- 
erned by either proliferation of domain walls [l8| or vor- 
tices UJ. Which one is realized, depends on the relative 
importance of the individual order parameter stiffness 
and the chiral locking term in the free energy. In the weak 
coupling Hubbard model it is likely that the chiral term 
is small, and thus the domain walls have lower energy 
cost. Then, one may expect that the finite-temperature 
transition is in the generalized Ising class 01 . 

In conclusion, we studied a model of itinerant elec- 
trons interacting with a chiral magnetic texture on a tri- 
angular lattice. Wc found that at 3/4 filling the system 
can become a spontaneous quantum Hall insulator, even 
without static magnetic ordering. One possible material 
realization is Nao.5Co02 [20| . 
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